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Abstract 

Major  impediments  to  the  wide-scale  implementation  of  hydrogen/air  fuel  cell  vehicles  are  the  lack  of  hydrogen  infrastructure  and  on-board 
hydrogen  storage.  One  proposed  source  of  hydrogen  exists  in  the  development  of  on-board  methanol  (and  other  hydrocarbon)  fuel  processors. 
Packaging  limitations  and  fuel  processor  performance  constraints  on  efficiency  and  transient  response  play  key  roles  in  vehicular  applications. 
These  constraints  may  be  addressed  by  considering  proper  thermal  integration  between  two  major  components  of  the  fuel  processor:  the 
reformer  and  catalytic  burner. 

The  focus  of  this  research  is  on  the  effects  of  the  catalytic  burner  on  reformer  performance  in  a  thermally  well-integrated  configuration. 
Specifically,  the  work  has  focused  on  the  generation  of  a  detailed  numerical  model  incorporating  kinetics  and  mass  and  heat  transfer  to 
accurately  characterize  the  burner.  Unlike  a  simple,  thermodynamic  model,  the  detailed  model  provides  a  level  of  complexity  necessary  to 
understand  the  impact  of  thermal  integration  on  reformer  transient  response,  reformate  composition,  and  emissions. 
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1.  Introduction 

Hydrogen/air  fuel  cell  vehicles  are  currently  receiving 
much  attention  as  an  alternative  to  internal  combustion 
engine  vehicles  (ICE  Vs)  because  of  their  potential  to  reduce 
fuel  consumption  and  emissions.  However,  the  combined 
problems  of  on-board  hydrogen  storage  and  the  lack  of 
hydrogen  infrastructure  represent  major  impediments  to 
their  wide-scale  adoption  as  replacements  for  ICE  Vs.  On¬ 
board  fuel  processors  that  generate  hydrogen  from  on-board 
liquid  methanol  (and  other  hydrocarbons)  have  been  pro¬ 
posed  as  alternative  hydrogen  sources  for  the  fuel  cell.  The 
indirect  methanol  fuel  cell  vehicle  fuel  processor  uses  steam 
reformation  of  methanol  to  produce  hydrogen  required  by 
the  fuel  cell.  Since  steam  reformation  is  an  endothermic 
process,  a  heat  source,  such  as  a  catalytic  burner,  is  required 
to  supply  the  necessary  thermal  energy. 

Packaging  limitations  and  fuel  processor  performance 
constraints  such  as  high  efficiency  and  fast  transient 
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response  play  key  roles  in  vehicular  applications.  Compo¬ 
nents  within  the  fuel  cell  system  have  transient  response 
times  on  the  order  of  milliseconds;  however,  the  fuel  pro¬ 
cessor  can  take  on  the  order  of  several  seconds  for  transient 
response.  These  constraints  may  be  addressed  by  considering 
proper  thermal  integration  between  the  reformer  and  burner. 

Several  configurations  of  the  burner  and  reformer  have 
been  proposed  and  have  been  published  in  the  literature 
[1-3]  such  as  shell-and-tube  heat  exchanger  and  concentric 
burner  and  reformer  arrangements.  However,  these  config¬ 
urations  have  limitations  in  terms  of  transient  response. 
Results  of  one  modeling  study  [1]  of  a  dynamic  reformer 
sized  to  produce  hydrogen  for  an  80  kWe  PEM  fuel  cell 
stack  show  ramp  down  response  greater  than  10  s  to  reach 
steady  state.  The  configuration  considered  is  a  series  of  tubes 
filled  with  copper  and  zinc  oxide  catalyst  surrounded  by  tube 
walls  and  hot  gas.  Similarly,  in  another  study  [2],  a  tubular 
reactor  was  modeled  that  generates  hydrogen  required  for  a 
400  kWe  PEM  fuel  cell  stack.  The  heat  is  transferred  to  the 
inner  wall  of  the  reformer  from  a  concentric  heat  source. 
Results  suggest  that  response  times  on  the  order  of  15-20  s 
are  possible.  Finally,  an  experiment  was  conducted  on  a  fuel 
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Nomenclature 

Unless  indicated  otherwise,  the  following  nomenclature 

applies  to  the  equations  presented  in  this  paper  Note  that 

“catalyst”  implies  catalyst  plus  the  catalyst  support  system. 

A 

reaction  surface  area  in  one  CSTR, 
(m  ;  product  of  Scat,  Ac,  £,  Ax) 

Ac 

reactor  cross-sectional  area  (m  ) 

Cor  [] 

Q 

concentration  (mol/m  ) 

Cp 

specific  heat  (J/mol  K) 

^Pcat 

catalyst  specific  heat  (J/kg  K) 

r 

^Pwall 

wall  specific  heat  (J/g  K) 

F  or  dn/dt  or  h 

flowrate  (mol/s) 

h 

enthalpy  (J/mol) 

hr> 

mass  transfer  coefficient  (m/s) 

ht 

heat  transfer  coefficient 
(W/m2  K) 

S 

open  frontal  area  fraction 

5cat 

l-S 

Methanol  reaction 

k\ 

gas  reaction  rate  constant  (s  ) 

k2 

surface  reaction  rate  constant  (m/s) 

Nitric  oxide  (NO)  reaction 

ku  k2 

reaction  rate  constants 

o 

(m  /(mol  s)) 

Pleat 

catalyst  mass  (kg) 

Plwall 

wall  mass  (g) 

n 

specie  quantity  (mol) 

r 

reaction  rate  (mol/s) 

R 

universal  gas  constant 
(8.315  J/mol  K) 

p 

pressure  (Pa) 

^reformer 

reformer  heat  rate 
requirement  (J/s) 

T 

temperature  (K) 

UA 

product  of  heat  transfer 
coefficient  and  area  (W/K) 

^burner— wall— side4wall 

product  of  the  heat  transfer 
coefficient  between  the  burner 
catalyst  and  wall  and  wall  area 

^reformer— wall— side^wall 

product  of  the  heat  transfer 
coefficient  between  the  wall  and 
reformer  catalyst  and  wall  area 

kyoid 

void  volume  of  one  CSTR, 

(m  ;  product  of  S ,  Ac,  Ax) 

Greek  symbols 

Ax 

length  of  one  CSTR  (m) 

A 

thermal  conductivity  (W/m  K) 

Pea 

catalyst  density  (kg/m  ) 

£ 

specific  surface  area 

(lrr/m  bed  volume 

(bed  volume  =  catalyst  +  void)) 

Subscripts 

i 

ith  specie 

g 

gas 

cat 

catalyst 

T 

total 

in 

entrance 

processor  system  in  which  catalytic  burners  were  placed  in 
an  annular  arrangement  at  the  bottom  of  a  50  kW  (H2) 
compact  reformer  [3].  Transient  response  times  were  mea¬ 
sured  to  be  about  20  s. 

Some  organizations  are  working  toward  improving  tran¬ 
sient  response  times  by  experimenting  with  and  modeling 
the  plate  configuration  [4-6],  typically  in  which  one  side  of  a 
plate  is  covered  with  reformer  catalyst  and  the  other  side  is 
covered  with  combustion  catalyst.  All  of  these  groups 
studied  units  that  could  produce  hydrogen  for  a  20  kWe 
fuel  cell  stack.  One  study  [4]  integrated  reforming,  combus¬ 
tion  and  selective  oxidation  (used  to  oxidize  the  CO  in  the 
reformate  which  can  otherwise  poison  the  fuel  cell)  catalysts 
coated  onto  aluminum  substrates  in  a  prototype  unit.  The 
heat  transfer  characteristics  of  the  compact  heat  exchange 
units  could  satisfy  transient  response  requirements  for  which 
the  goal  was  set  at  less  than  5  s. 

The  literature  presents  examples  of  on-going  experimen¬ 
tal  and  modeling  efforts  in  characterizing  a  thermally  well- 
integrated  system.  While  studies  have  been  conducted  on  the 
effects  of  a  thermally  well-integrated  system  on  reformer 
performance  under  steady  state  conditions,  this  paper  pro¬ 
vides  simulation  results  of  a  detailed  model  run  under 
varying  and  dynamic  conditions.  The  model  was  developed 
using  Matlab/Simulink  software.  The  results  of  research 
discussed  in  this  paper  are  based  on  this  model,  which  is 
part  of  a  fuel  cell  system  designed  for  integration  into  a 
larger  indirect  methanol  fuel  cell  vehicle  model. 

This  paper  is  comprised  of  two  major  sections:  (1)  a 
description  of  the  steps  taken  in  generating  a  detailed  model 
that  will  accurately  characterize  a  thermally  integrated 
burner/reformer  system,  and  (2)  a  description  of  simulation 
results  that  show  the  effects  of  thermal  integration  on:  (a) 
reformer  temperature  and  emissions  especially  when  vary¬ 
ing  spatial  parameters  such  as  catalyst  loading,  (b)  sensi¬ 
tivity  to  changes  in  input  flow  or  thermal  mass,  (c)  dynamic 
response,  and  (d)  burner  emissions.  Vehicle  startup  and 
control  issues  and  discussion  of  a  detailed  reformer  model 
are  not  part  of  this  paper.  A  paper  on  the  detailed  reformer 
model  is  forthcoming. 

2.  Background  of  the  detailed  model 

2.7.  Indirect  methanol  fuel  cell  (IMFC)  system 

The  burner  is  part  of  the  fuel  processor  that  is  part  of  a 
larger  fuel  cell  system.  As  shown  in  Fig.  1,  the  burner 
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Fig.  1.  Schematic  of  the  IMFC  system  as  studied  in  the  IMFC  vehicle  model. 


receives  methanol  from  an  on-board  tank,  air  from  a  blower 
and,  in  some  operating  conditions,  exhaust  from  the  fuel  cell 
anode,  which  consists  of  hydrogen,  C02  and  water  vapor. 
The  burner  provides  heat  to  a  preheater  and  the  steam 
reformer. 

Typically,  the  burner,  via  an  intermediate  heat  transfer 
medium,  heats  the  reformer.  In  the  simplest  shell  and  tube 
model,  the  burner  hot  gases  provide  the  energy  needed  to 
heat  the  reformer.  As  high  temperatures  are  not  desired  due 
to  NO*  formation  and  gases  have  low  heat  capacity,  the 
flowrate  of  the  gases  must  be  very  high  to  ensure  adequate 
heat  transfer  to  the  reformer.  High  temperatures  and  flow- 
rates  are  not  desirable  operating  conditions,  and  alternate 
methods  such  as  hot  oil  as  an  intermediate  thermal  medium 
have  been  used  in  some  cases.  This  method  has  its  own 
drawbacks  since  hot  oil  would  introduce  additional  thermal 
mass  that  would  need  to  be  heated.  Due  to  the  limitations  of 
conventional  approaches,  a  survey  of  existing  on-going 


Fig.  2.  Physical  representation  of  thermally  integrated  burner  and  reformer 
bi-catalyst  plate  configuration. 


industrial  efforts  in  this  area  indicates  that  well-integrated 
burner/reformer  designs  are  currently  being  pursued  by  a 
number  of  different  organizations  [4-6].  One  abstracted 
physical  representation  of  such  a  model  is  as  shown  in  Fig.  2. 

One  side  of  each  plate  is  coated  with  the  reformer  catalyst 
and  the  other  side  with  the  burner  catalyst.  This  burner/ 
reformer  configuration  uses  conduction  across  the  width  of  a 
thin  plate  as  a  mode  of  heat  transfer  from  the  burner  to  the 
reformer.  This  results  in  very  rapid  transfer  of  heat  from  the 
burner  to  the  reformer  compared  to  some  alternate  methods. 
Even  though  improvements  in  other  burner/reformer  con¬ 
figurations  continue,  the  plate  design  is  chosen  as  the 
physical  burner/reformer  model  in  this  paper. 

2.2.  Mathematical  modeling 

In  the  process  of  developing  the  model,  the  level  of  detail 
had  to  be  assessed.  A  simple,  thermodynamic  model  with 
inlet  and  exit  states  would  provide  only  such  information  as 
adiabatic  temperature  and  exit  composition  at  equilibrium. 
However,  there  are  some  reasons  why  a  more  detailed  model 
is  required.  First,  within  the  fuel  cell  system  the  fuel 
processor  has  a  slow  transient  response  relative  to  the  other 
components  and  it  cannot  be  assumed  that  thermal  equili¬ 
brium  is  reached  for  every  input  flow  change.  Only  a 
chemical  kinetic  model  can  determine  the  instantaneous 
output  conditions  for  a  dynamic  input.  Second,  as  the  burner 
and  reformer  are  thermally  integrated,  the  instantaneous 
burner  conditions  directly  affect  the  reformer  conditions. 
Finally,  the  literature  provided  sources  of  kinetic  data  which 
could  be  used  to  validate  the  model  against  laboratory 
results.  Developing  a  detailed  model  allowed  a  level  of 
accuracy  to  be  ascertained. 

Once  it  was  determined  that  a  kinetic  model  should  be 
used,  the  modeling  technique  had  to  be  established.  For 
the  catalytic  burner,  a  plug  flow  reactor  configuration 
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Fig.  3.  Ten  CSTR  discretization  of  plate  configuration. 


represented  by  a  series  of  continuously  stirred  tank  reac¬ 
tors  (CSTRs)  was  assumed.  Ideally,  a  plug  flow  reactor  has 
no  mixing  of  the  flow  in  the  axial  direction  and  perfect 
mixing  in  the  radial  direction.  Concentration  and  tempera¬ 
ture  gradients  occur  in  the  axial  direction  only.  Each 
CSTR,  however,  has  perfect  mixing  such  that  the  reactor 
contents  are  spatially  uniform.  The  modeling  methodology 
is  to  discretize  the  reactor  into  several  elements  in  the 
axial  direction,  and  assume  all  the  quantities  to  be  uniform 
in  each  element  in  a  finite  difference  approximation  of 
spatial  gradients.  The  number  of  elements  is  dictated  by 
the  balance  between  accuracy  and  computational  com¬ 
plexity.  Since  the  reformer  model  was  also  constructed 
using  the  CSTR  approach  and  the  burner  was  integrated 
with  the  reformer,  it  was  determined  that  10  CSTRs  would 
be  used  based  on  a  literature  search  conducted  for  the 
reformer  [2,7].  To  relate  the  plate  configuration  to  the  10 
CSTR  modeling  methodology,  Fig.  3  is  shown  with  lines 
to  illustrate  discretization. 

2.3.  Burner  kinetics 


substrate 


catalyst  and  wash  coat 


open  frontal  area 


monolith  channel 


Fig.  4.  Schematic  representation  of  a  monolith  channel  (not  to  scale). 


The  reactor  is  assumed  to  operate  at  atmospheric  pressure, 
the  gas  mixture  is  an  ideal  gas,  and  pressure  drop  is 
neglected. 

2.3.1.  Rate  expressions  used 

Rate  expressions  are  provided  in  [8]  for  the  gas  and 
surface  reaction,  where  the  reaction  rate  constants  (k{  and 
k2 ,  described  below)  were  experimentally  determined  and 
one-step  Arrhenius  rate  expressions  for  complete  oxidation 
reactions  for  the  surface  and  gas  were  developed.  Based  on 
the  earlier  work  of  [8]  of  methanol  oxidation  on  platinum,  it 
was  determined  that  the  reaction  order  of  methanol  was  first- 
order  and  oxygen  was  zero-order.  The  methanol  reaction 
rates  were  determined  under  steady  state  conditions  by  the 
difference  between  the  methanol  concentration  at  the  inlet 
and  the  outlet  of  the  experimental  catalytic  reactor.  The 
methanol  oxidation  equation  is  shown  in  Eq.  (1). 

CH3OH  + §02  =  C02  +  2H20  (1) 

The  homogeneous  reaction  rate  of  methanol  is  shown  in 
Eq.  (2). 

rCH3OH  =  CcH3OHg  ^void  (2) 

CcH3OHg  is  the  concentration  of  methanol  in  the  gas,  reaction 
rate  constant  k\  =  e*-6-1*10  /(burner) +3. 2)^  anc[  yvoid  the 

burner  void  volume  for  one  CSTR,  which  is  calculated  using 
the  open  frontal  area. 

As  shown  in  Eq.  (3)  the  concentration  of  surface  metha¬ 
nol,  CcH3OHcat,  had  to  be  determined  by  equating  two  surface 
reaction  rate  expressions: 

^CH3OH  =  (CcH3OHg  —  CcH3OHcatM> 

rCH3OH  =  £2CcH3OHcat7l  (3) 


The  plug  flow  configuration  for  the  burner  was  modeled 
based  on  the  experimental  and  numerical  work  done  by  [8]. 
The  burner  kinetics  are  based  on  a  platinum-loaded  gamma 
alumina  monolith  reactor  and  assume  oxidation  of  methanol 
in  a  one- step  reaction.  Fig.  4  illustrates  a  representation  of  a 
monolith  channel  where  reactions  take  place  in  the  gas  and 
on  the  catalyst  particles  which  are  washcoated  onto  the 
substrate  surface.  The  open  frontal  area,  which  is  the 
porosity,  void,  or  free  cross-sectional  area  available  for  flow, 
is  64%.  It  was  determined  in  the  present  analysis  that  a  total 

Q 

reactor  volume  of  0.0041  m  is  required  to  ensure  complete 
conversion  of  the  methanol. 


Solving  for  CcH3OHcat  yields  the  following  expression: 


CCH3OHcat 


6^CH3OHg^D 
k2  +  hj) 


The  term  /zD  is  the  mass  transfer  coefficient,  the  reaction  rate 
constant  k2  =  e(-5-9xlo4/(^rbumer)+22.45)  when  rburner  <  380  K 
and  k2  =  e^-64  xl03/(RTbumer)+5.8 )  wJien  380  K  <  Tburner  < 
830  K,  and  A  is  the  reaction  surface  area.  The  reaction  rates 
of  the  other  species  in  the  methanol  reaction  are  based  on 
Eqs.  (2)  and  (3)  and  the  reaction  stoichiometry.  The  mass 
transfer  coefficient,  /zD,  is  0.235  m/s  and  was  determined 
from  [9]. 


M.  Sundaresan  et  al.  /  Journal  of  Power  Sources  113  (2003)  19-36 


23 


Energy  and  species  balance  calculations  are  performed  in 
each  CSTR.  Eqs.  (4)— (8)  are  used  in  the  species  balance. 


d  t 


The  term  d «/cat  / d t  is  the  flowrate  due  to  reactions  and  ricat  the 
reaction  rate,  both  for  the  ith  surface  specie.  Note  the  sign 
convention  for  reaction  rates  as  used  in  the  model. 


—jf  =  ^i’(in)  —  (rig  +  ricat)  —  Fi  (5) 

The  term  dnig  / d t  is  the  flowrate  due  to  reactions  and  flows 
into  and  out  of  the  burner,  rig  the  reaction  rate,  the 
flowrate  into  the  burner  and  Fi  the  flowrate  out  of  the  burner, 
all  for  the  ith  gas  specie. 


Tlj 


PV void 


FI burner 


The  term  nT  is  the  total  moles,  assuming  ideal  gas  condi¬ 
tions,  P  the  burner  total  pressure,  Tbumer  the  burner  tem¬ 
perature,  and  R  the  universal  gas  constant.  Note  that 
calculations  are  done  with  moles  and  therefore  use  the 
universal  gas  constant  R ;  however,  while  moles  are  not 
conserved,  mass  is  conserved.  This  condition  was  verified 
within  the  model.  The  initial  conditions  for  moles  are 
no 2  =  0.2 nj,  and  n^2  =  0.8 nT  assuming  only  air  exists 
within  the  burner  at  t  =  0. 


The  term  nt  is  the  moles  for  the  ith  specie  and  FT  is  the  total 
flowrate  which  is  further  defined  in  Eq.  (8).  The  initial 
condition  for  Tbumer  is  550  K. 


Ft  =  5>0n)  -  yyng  +  rij  + 


P y void  d 7bumer 


dt2 

/v  1  burner 


d  t 


(8) 


The  terms  included  in  the  energy  balance  equations  depend 
on  the  level  of  required  complexity  to  accurately  and 
realistically  model  the  burner.  Two  plug  flow  models  were 
examined:  (1)  a  heterogeneous  model,  which  includes  heat 
transfer  between  the  surface  and  gas  and  conduction  in  the 
substrate  along  the  burner,  i.e.  surface  and  gas  temperatures 
are  different  and  (2)  a  pseudo-homogeneous  model,  which 
assumes  that  the  gas  and  surface  are  at  the  same  temperature. 
In  both  cases,  radiation  within  the  combustor  and  to  the 
surroundings  is  neglected.  In  the  thermally  integrated  case, 
temperatures  are  lower  than  1000  K.  In  combustors  where 
temperatures  are  above  1000  K  and  can  reach  1500  K,  [10] 
recommends  that  heat  transfer  by  radiation  be  considered. 
Heat  loss  from  the  combustor  is  assumed  only  to  the 
reformer  and  not  to  the  ambient  environment. 

Simulations  were  run  to  compare  the  results  of  the 
heterogeneous  and  pseudo-homogeneous  model  to  deter¬ 
mine  the  importance  of  these  energy  balance  differences. 
In  a  full  vehicle  model,  detailed  components  models  such 
as  for  the  heterogeneous  burner  would  require  hours  of 


computation  time,  so  the  comparison  was  necessary  to 
determine  whether  a  simpler  model  such  as  the  pseudo- 
homogeneous  model  could  be  used.  The  results  of  the 
comparison  and  the  energy  and  species  balance  equations 
are  provided  in  the  next  section. 

2.4.  Thermal  coupling 

The  method  used  to  thermally  couple  the  burner  and 
reformer  CSTRs  lies  in  linking  the  energy  balance  equa¬ 
tions.  Since  heat  transfers  from  the  burner  catalyst  to  the  wall 
and  the  wall  to  the  reformer  catalyst,  energy  balance  equa¬ 
tions  were  required  for  the  three  sections.  However,  it  was 
necessary  to  determine  the  importance  of  heat  transfer 
between  the  surface  and  gas  and  conduction  in  the  substrate 
along  the  burner,  which  are  modes  present  in  the  hetero¬ 
geneous  burner  configuration.  Such  a  model  would  account 
for  a  difference  in  surface  and  gas  temperatures;  however, 
the  pseudo-homogeneous  model  assumes  the  surface  and 
gas  temperatures  are  the  same.  As  shown  in  the  following 
figures,  the  comparison  between  the  heterogeneous  and 
pseudo-homogeneous  models  reveals  little  difference,  espe¬ 
cially  in  the  reformer  exit  hydrogen,  which  is  the  key 
parameter  in  question.  Note  that  the  results  of  validation 
runs  for  the  heterogeneous  model  compared  to  results  of  [8] 
are  shown  in  Appendix  A.  Also  note  that  a  paper  on  the 
detailed  reformer  model  is  forthcoming. 

Fig.  5  shows  a  comparison  of  reformer  exit  hydrogen  for 
the  two  burner  models. 

The  simulations  were  run  with  a  methanol  input  flowrate 
stepped  up  to  0.35  mol/s  for  the  reformer  (and  input  steam 
at  a  steam/methanol  ratio  of  1.3)  and  0.05  mol/s  for  the 
burner  (and  input  air  at  a  relative  air/fuel  ratio  of  1.5).  These 
flowrates  represent  one  particular  operating  condition  for  the 
fuel  processor.  For  reference,  the  reaction  for  the  endother¬ 
mic  steam  reformation  of  methanol  is  shown  in  Eq.  (9). 
Results  show  only  1  %  difference  in  reformer  exit  hydrogen, 
which  is  considered  negligible. 

CH3OH  +  H20  =  C02  +  3H2  (9) 


1500  16(H)  17(H)  1800  19(H)  20(H) 

time,  seconds 


Fig.  5.  Comparison  of  reformer  exit  hydrogen.  ( - )  Pseudo- 

homogeneous  model;  ( — )  heterogeneous  model.  Note  scaling;  see  Fig.  6 
for  full  step  response. 
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Fig.  6.  Step  response  of  exit  hydrogen.  ( - )  Pseudo-homogeneous 

model;  ( — )  heterogeneous  model. 


An  explanation  of  the  transient  response  follows.  In  Fig.  6 
the  step  occurs  at  1500  s  and  there  is  the  appearance  of  an 
instantaneous  response. 

The  reason  for  this  behavior  is  because  the  reformer  has  a 
finite  volume  and  the  latter  half  of  the  fuel  processor  is  rich 
in  hydrogen.  When  the  step  input  of  methanol/water  mixture 
comes  into  the  reformer,  it  effectively  pushes  the  hydrogen- 
rich  reformate  out  of  the  reformer  (note  that  the  reformer  is 
operated  in  a  constant  pressure  mode).  However,  this  sce¬ 
nario  cannot  be  sustained  and  the  response  follows  first- 
order  behavior  until  it  approaches  steady  state. 

Based  on  the  comparison  results,  it  was  decided  to  use  the 
simpler,  pseudo-homogeneous  energy  balance  as  shown  in 
Eq.  (10).  The  main  parameter  used  in  the  models  to  thermally 
couple  the  burner  and  reformer  is  UA ,  where  U  is  the  heat 
transfer  coefficient  and  A  is  the  area.  Specifically,  UA  is  the 
product  of  the  heat  transfer  coefficient  between  the  burner 
catalyst  and  wall  and  wall  area  (f/bumer-waii-side^waii)  and  the 


product  of  the  heat  transfer  coefficient  between  the  wall  and 
reformer  catalyst  and  wall  area  (t/reformer-waii-sideAwaii).  The 
“thermal  coupling”  is  made  as  follows:  rbumer,  or  the  burner 
temperature  in  Eq.  (10)  is  calculated  and  used  in  the  wall 
energy  balance  in  Eq.  (11).  rwan,  or  the  wall  temperature  in 
Eq.  (11)  is  calculated  and  used  in  the  heat  transfer  term  of 
Eq.  (12).  The  heat  transfer  term  in  Eq.  (12)  equals  the  heat 
the  burner  must  provide  and  is  used  for  the  term  <2reformer  in 
Eq.  (10).  Note  that  “catalyst”  implies  catalyst  plus  the 
catalyst  support  system. 


^stored 


\  d7bumer 


^  ^  T)(in)^(in)  ^reformer  ^  ^  T 


drii 

dt 


(10) 


mwall  CPwa|| 


dT, 


wall 


dt 


U burner — wall — side  A  wall  (  Tburner  -  Twan) 


(ID 


dTref 


(mcat  Cpcat  +  ^  rii  Cp . )  ^ 

^reformer— wall— sideA wall  (Twall  Tref) 


(12) 


The  initial  condition  for  the  temperatures  are  Tburner  = 
TWaii  =  Tref  =  550  K,  where  the  burner,  wall  and  reformer 
are  assumed  to  be  fully  warmed  up  at  time  t  =  0.  A 
representation  of  the  thermal  coupling  between  the  10 
CSTRs  of  the  reformer  and  burner  is  illustrated  in  Fig.  7. 

f^burner— wall— sideAwall  dnd  ^7reformer— wall— sideA wall  had  tO  be 


determined.  However,  it  soon  became  apparent  that  the  fuel 
processor  would  need  to  be  scalable  using  the  UA  para¬ 
meters  while  maintaining  the  validity  of  using  the  existing 
kinetic  data. 


Fig.  7.  Illustration  of  thermal  coupling. 
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2.5.  Scaling 

Since  kinetic  expressions  are  non-linear,  scaling  a  tubular 
reactor  is  not  a  linear  process.  For  example,  if  an  increase  in 
output  flowrate  of  X  times  is  required,  increasing  reactor  size 
either  by  catalyst  mass  or  volume  by  X  times  would  decrease 
the  time  it  takes  to  get  to  the  final  product,  but  not  increase 
the  amount  of  the  final  product.  Reaction  rates  are  typically 
per  unit  catalyst  mass  or  volume,  so  increasing  either 
catalyst  mass  or  volume  simply  increases  the  rate  at  which 
a  reactant  is  consumed  or  product  is  formed.  Fig.  8  illustrates 
three  plots  of  response  time. 

Plot  (a)  represents  the  response  of  a  25  kg  reformer,  used 
as  a  baseline.  Plot  (bl)  represents  the  response  of  a  35  kg 
reformer.  Increasing  the  catalyst  mass  allows  the  reaction  to 
proceed  more  quickly  to  the  same  steady  state  exit  flow. 
However,  as  seen  in  plot  (b2),  which  also  represents  the 
response  of  a  35  kg  reformer,  increased  catalyst  mass  can 
take  more  time  to  reach  steady  state  if  thermal  energy  is 
inadequate. 

Increasing  the  flowrate  by  X  times,  without  an  increase  in 
catalyst  mass  or  volume,  would  result  in  lower  conversion. 
Fig.  9  shows  the  effect  of  increasing  input  flow  on  the  same 
reactor.  Plots  (a)  and  (b)  represent  flowrates  of  0.05  and 
0.5  mol/s,  respectively.  The  exit  conversion  is  lower  for  the 
higher  input  flowrate. 

The  plate  configuration  is  effectively  having  several 
reactors  in  parallel;  therefore,  scaling  is  possible  by  increas¬ 
ing  the  number  of  plates  as  shown  in  Fig.  10.  Based  on  the 
maximum  required  exit  hydrogen,  one  could  manifold  an 
input  flowstream  to  the  correct  number  of  plate  channels 
such  that  individual  channel  flows  do  not  exceed  the  max¬ 
imum  reactant  flow  that  still  achieves  the  desired  conversion 
given  the  reaction  kinetics.  It  is  assumed  that  the  plates  are 
fully  warmed  up.  Furthermore,  it  is  assumed  that  all  plates 
are  in  use  during  system  operation.  Even  if  some  plates  are 
warm  but  not  in  use,  i.e.  no  flow  in  the  plate  channels 


Fig.  8.  Effect  of  increasing  catalyst  mass  on  response  time:  (a)  25  kg, 
baseline,  (bl)  35  kg,  faster  thermal  response  than  baseline,  (b2)  35  kg, 
slower  thermal  response  than  baseline. 


0  100  200  300  400  500  600 


time,  seconds 


Fig.  9.  Effect  of  flowrate  on  conversion:  (a)  0.05  mol/s,  >99%  conversion, 
(b)  0.5  mol/s,  ~78%  conversion. 


because  full  power  operation  is  not  yet  needed,  the  response 
may  be  different  than  if  all  plate  channels  had  a  uniform, 
lower  flowrate  before  being  brought  up  to  full  power. 

While  the  purpose  of  constructing  and  using  this  model  is 
not  to  design  systems,  as  this  is  dependent  on  vehicle  design, 
it  is  important  to  be  realistic  in  model  parameters  and  make  a 
“best-guess”  assessment  of  packaging  considerations.  To 
make  scaling  in  the  model  parameters  flexible  and  realistic 
for  the  plate  configuration,  it  was  required  that  the  UA 
parameters  become  independent  of  individual  plate  area. 
The  best  way  to  determine  these  UA  values  would  be  from 
microscale  experiments.  It  might  also  be  possible  to  estimate 
UA  from  geometrical  and  thermal  properties  or  alternatively, 
given  a  value  for  UA ,  the  correct  set  of  thermal  properties 
can  lead  to  initial  geometry  estimates. 

In  the  initial  part  of  the  analysis,  the  values  for  UA  were 
selected  using  the  work  done  earlier  by  [7] .  That  study  on  a 
shell  and  tube  reformer/burner  configuration  illustrated  the 
need  for  substantially  higher  heat  transfer  coefficients  to 
achieve  exit  H2  response,  from  10  to  90%,  of  the  order  of  2  s, 
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Fig.  10.  Scalability  with  plate  configuration. 
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which  is  the  response  time  required  by  the  fuel  cell  system 
specifications  used  [11].  Therefore,  the  initial  studies  of  the 
thermally  integrated  reformer/burner  system  were  based  on 
an  assumption  of  setting  both  t/bumer-waii-sidAvaii  and  Ure. 
former— wall— sideAwaii  equal  to  60  W/K  per  plate  (for  100  plates), 
which  is  an  order  of  magnitude  higher  than  the  values 
studied  by  [7].  A  more  detailed  discussion  of  the  UA 
parameters  can  be  found  in  Appendix  B. 

3.  Model  studies  and  results — thermal 
integration  effects 

As  mentioned  earlier,  some  advantages  of  the  discretiza¬ 
tion  of  the  reformer  and  burner  into  ten  elements  include  the 
ability  to  see  the  effects  of  parametric  changes  and  obtain 
local  values  at  every  point  along  the  reactor.  Discretization 
allows  for  thermally  integrating  the  burner  and  reformer  and 
varying  parameters  to  observe  the  interaction  between  the 
two  components.  For  example,  variations  have  been  studied 
for  the  burner  catalyst  loading  and  wall  thermal  mass.  A 
four-section  discussion  describes  simulation  results  that 
show  the  effects  of  thermal  integration  on:  (a)  reformer 
temperature  and  emissions  especially  when  varying  spatial 
parameters  such  as  catalyst  loading  (Section  3.1),  (b)  sensi¬ 
tivity  to  changes  in  input  flow  or  thermal  mass  (Section  3.2), 
(c)  dynamic  response  (Section  3.3),  and  (d)  burner  emissions 
(Section  3.4).  Note  that  despite  the  model’s  ability  to  accept 
anode  exhaust  as  shown  in  Fig.  1,  simulation  results  in  this 
section  are  based  on  methanol  input  alone  to  observe  the 
effects  of  thermal  integration  and  eliminate  the  complexity 
of  an  additional  fuel. 

3.1.  Catalyst  loading 

Varying  the  burner  catalyst  loading  plays  an  important 
role  in  reducing  reformer  steady  state  temperature  and 
carbon  monoxide  output.  Reducing  reformer  temperature 
protects  the  catalyst  from  operating  in  a  region  in  which  it 
can  deactivate  due  to  sintering.  The  reformer  CO  formation 
is  a  strong  function  of  temperature  based  on  the  reformer 
kinetic  model  used  [12].  Reducing  the  reformer  temperature 
corresponds  to  a  reduced  amount  of  CO. 

Carbon  monoxide  can  be  a  significant  problem  from  the 
viewpoint  of  the  fuel  cell  stack.  Even  while  there  is  a  CO 
cleanup  system  in  place,  minimizing  or  eliminating  CO  from 
the  reformer  would  be  an  ideal  condition.  Modifying  para¬ 
meters  within  the  reformer  to  control  or  optimize  CO  output 
would  need  to  be  considered,  and  as  described  below,  these 
modifications  can  yield  favorable  results. 

3.1.1.  Simulation  setup 

Close  thermal  integration  results  in  different  responses  to 
different  types  of  burner  catalyst  loading.  Three  configura¬ 
tions  along  the  burner  are  discussed  and  are  shown  in  Fig.  1 1 : 
(a)  uniform  loading,  (b)  gradually  increasing  loading,  and 


Fig.  11.  Three  catalyst  loading  distribution  schemes  along  burner:  (a) 
uniform,  (b)  gradually  increasing,  (c)  high  concentration — second  half  of 
burner. 

(c)  high  concentration  of  catalyst  in  the  second  half  of  the 
burner.  In  all  cases  there  is  the  same  total  amount  of  catalyst. 
Note  that  the  reformer  catalyst  loading  is  uniform  in  this 
analysis.  The  following  parameters  were  analyzed:  (1) 
burner  temperature,  (2)  reformer  temperature,  (3)  reformer 
CO  flow,  and  (4)  reformer  output  H2. 

3.1.2.  Simulation  results 

Figs.  12  and  13  show  the  effects  of  the  three  catalyst 
loading  scenarios  on  the  burner  and  reformer  temperature 
profiles.  Configuration  (a)  represents  the  uniform  loading 
scenario  which  has  the  highest  temperature  profile.  Config¬ 
uration  (b)  is  the  gradually  increasing  loading  scenario  for 
which  the  temperature  in  the  initial  CSTRs  is  slightly  lower. 
For  configuration  (c),  the  temperature  hovers  around  550  K 
with  the  exception  of  the  jump  from  the  fifth  to  the  sixth 
CSTR  at  which  point  the  amount  of  burner  catalyst 
increases,  increasing  the  local  reaction  temperature. 

The  reformer  CO  formation  is  a  strong  function  of 
temperature  and  this  relationship  is  reflected  in  Fig.  14 
which  shows  lower  steady  state  molar  flowrates  of  CO  for 
configurations  (b)  and  (c). 


Fig.  12.  Burner  temperature  for  varied  catalyst  loading:  (a)  uniform,  (b) 
gradually  increasing,  (c)  high  concentration — second  half  of  burner. 
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Fig.  13.  Reformer  temperature  for  varied  catalyst  loading:  (a)  uniform,  (b) 
gradually  increasing,  (c)  high  concentration — second  half  of  burner. 


In  Fig.  14,  while  the  CO  flowrate  in  configuration  (c) 
appears  to  continue  to  increase  at  the  end  of  the  reformer, 
Fig.  15  shows  that  the  output  H2  ceases  to  change  at  the  end 
of  the  reformer  and  illustrates  how  the  same  output  hydrogen 
in  configurations  (a)  and  (b)  could  be  achieved  with  the 
lower  CO  penalty  of  configuration  (c). 

While  many  catalyst  distributions  may  not  have  an 
impact  or  may  even  have  an  adverse  impact  on  reformer 
performance,  the  results  for  the  example  loading  schemes 
used  in  this  analysis  show  that  varied  catalyst  loading  could 
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Fig.  14.  Reformer  CO  output  for  varied  catalyst  loading:  (a)  uniform,  (b) 
gradually  increasing,  c)  high  concentration — second  half  of  burner. 


play  a  role  in  decreasing  reformer  temperature  and  CO 
emissions. 

3.2.  Sensitivities 

Thermal  integration  makes  the  system  more  sensitive  to 
perturbations  in  the  input  flow  stream.  The  current  model 
assumes  no  time  lag  between  the  methanol  sent  into  the 
reformer  and  the  methanol  sent  into  the  burner  to  generate 
the  necessary  heat  for  the  reformer  reactions  to  occur.  How¬ 
ever,  a  time  lag  between  these  flows  can  affect  the  temperature 
distribution  of  the  thermally  integrated  components,  possibly 
jeopardizing  reformer  catalyst  integrity,  decreasing  methanol 
conversion  and  increasing  carbon  monoxide  production.  In 
this  section,  different  scenarios  of  thermal  integration  are 
presented.  At  the  same  time,  the  impact  of  a  lag  on  a  thermally 
integrated  system  is  investigated. 

3.2.1.  Impact  of  thermal  integration  on  reformer 
and  impact  of  lag  on  thermal  integration 

Based  on  the  fuel  cell  stack  H2  requirement,  the  controller 
will  send  the  required  amount  of  methanol  to  the  reformer  and 
based  on  the  heat  requirement  of  the  reformer,  the  controller 
will  send  the  required  amount  of  methanol  to  the  burner. 
However,  some  drawbacks  exist,  including  possible  reformer 
catalyst  sintering  due  to  a  higher  reformer  temperature  as  a 
result  of  the  thermal  integration.  Furthermore,  the  increase  in 
reformer  temperature  accelerates  the  production  of  carbon 
monoxide  under  certain  transient  conditions  since  the  refor¬ 
mer  CO  formation  is  a  strong  function  of  temperature  based 
on  the  reformer  kinetic  model  used. 

Additionally,  this  situation  can  be  exacerbated  by  a  con¬ 
dition  in  which  the  methanol  flowing  to  the  reformer  leads 
(or  lags)  the  methanol  flowing  into  the  burner  to  generate  the 
required  thermal  energy  for  the  steam  reformation  to  occur. 
A  shift  could  cause  the  burner  temperature  to  be  higher  than 
that  required  by  the  reformer  in  a  particular  CSTR,  jeopar¬ 
dizing  reformer  catalyst  integrity,  causing  less  methanol  to 
be  converted  into  hydrogen  or  more  carbon  monoxide  to  be 
generated. 
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Fig.  15.  Reformer  H2  output  for  varied  catalyst  loading:  (a)  uniform,  (b) 
gradually  increasing,  (c)  high  concentration — second  half  of  burner. 


3.2.2.  Simulation  setup 

In  this  analysis,  an  arbitrarily  selected  2  s  lag  in  burner 
step  input  methanol  flow  was  investigated.  Preliminary 
results  for  a  lag  to  a  step  response  are  discussed,  and  show 
that  there  appears  to  be  minimal  effect  of  a  lag  of  2  s.  Several 
parameters  were  analyzed: 

•  Burner  temperature 

•  Reformer  temperature 

•  Reformer  wall  temperature 

•  Reformer  MeOH 

•  Reformer  CO 

•  Output  H2 

The  relative  differences  between  the  lag  case  and  ideal 
case  are  small  but  more  pronounced  for  the  reformer  MeOH, 
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Fig.  16.  Step  input  flowrate  for  sensitivity  analysis.  Reformer  ( — ),  burner 
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Fig.  18.  Reformer  exit  MeOH — Scenarios  1  and  2;  Scenario  2  is  a  more 
closely  thermally  integrated  system  than  Scenario  1.  Ideal  condition  ( — ), 
lagging  condition  (...). 


CO  and  H2  than  the  temperatures;  therefore,  only  plots  of  the 
flows  are  shown. 

In  the  step  input  analysis  the  flows  are  compared  using 
three  scenarios.  Scenario  1  is  the  baseline  for  which  Uburner_ 

wall— side^wall  dnd  Cj-eformer— wall— side^wall  dre  CCJUal  tO  60  W/K 

per  plate  (for  100  plates).  It  was  compared  with  two  other 
scenarios  in  which  the  thermal  coupling  was  enhanced.  For 
Scenario  2,  both  UA  values  are  10  times  the  baseline.  For 
Scenario  3,  the  UA  values  are  the  same  as  in  Scenario  2; 
however,  the  wall  mass  is  reduced  10-fold.  As  discussed  in 
the  step  input  section,  the  increased  thermal  coupling  can 
improve  overall  methanol  conversion.  However,  the  lower 
thermal  mass  in  Scenario  3  increases  the  sensitivity  to  the  lag 
in  burner  input  methanol. 

In  the  model,  the  delay  was  introduced  in  the  controller 
before  the  burner  methanol  input.  Simulations  were  run  with 
a  2  s  input  delay.  Fig.  16  shows  a  snapshot  in  time  of  the 
methanol  step  input  used  for  this  analysis.  The  step  input  is 
preceded  by  an  idle  condition  for  200  s.  The  reformer  and 
burner  flowrates  reach  a  maximum  of  0.5  and  0.05  mol/s, 
respectively. 

3.2.3.  Simulation  results 

The  reformer  exit  unconverted  methanol  is  shown  in 
Fig.  17  versus  time  in  seconds  in  response  to  the  input 


time,  seconds 

Fig.  17.  Reformer  exit  MeOH — Scenario  1.  Ideal  condition  ( — ),  lagging 
condition  (...). 


shown  in  Fig.  16.  The  plot  represents  the  unconverted 
methanol  normalized  by  the  input  methanol.  There  is  a 
slight  increase  in  unconverted  methanol  in  the  lagging 
condition,  which  is  an  indication  that  fewer  reactions  are 
occurring  in  the  reformer  as  a  result  of  the  lag  in  burner  input 
methanol. 

Furthermore,  the  initial  spike  in  unconverted  methanol  for 
both  conditions  represents  the  maximum  amount  of  methanol 
in  the  reactor  at  the  time,  which  corresponds  to  the  step  up  to 
the  maximum  input  flowrate.  The  first-order  response  is  due 
to  the  reformer  thermal  mass,  which  maintains  a  temperature 
at  which  activity  can  occur  and  methanol  can  be  converted 
(note:  the  reformer  methanol  flow  would  not  be  a  step 
response  when  used  as  a  part  of  an  overall  system.  It  is  used 
more  as  a  diagnostic  tool  for  this  and  other  related  studies). 

An  interesting  point  can  be  illustrated  regarding  the  effect 
on  exit  methanol  in  a  thermally  integrated  system.  In  further 
improving  thermal  integration  by,  for  example,  increasing 

the  f^burner— wall— side^w all  dnd  ^reformer— wall— side^w all  1 0-fold, 

the  effects  of  closer  thermal  coupling  are  more  pronounced. 
Fig.  18  shows  the  effects  of  thermal  integration  on  uncon¬ 
verted  methanol,  which  decreases  for  the  more  closely 
thermally  integrated  system  (Scenario  2). 

Fig.  19  illustrates  exit  MeOH  for  Scenario  3  to  be 
compared  with  Scenarios  1  and  2  in  the  previous  figure. 


time,  seconds 


Fig.  19.  Reformer  exit  MeOH — Scenario  3.  Ideal  condition  ( — ),  lagging 
condition  (...). 
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Fig.  20.  Reformer  exit  temperature  before  and  after  step  input. 
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Fig.  22.  Reformer  exit  CO — ideal  condition. 


Scenario  2  represents  improved  methanol  conversion.  The 
exit  MeOH  after  the  stepup  for  Scenario  3  exceeds  that  of 
Scenarios  1  and  2  due  to  reduced  reaction  rates  resulting 
from  lower  exit  temperatures.  The  spike  in  exit  MeOH  for  all 
scenarios  corresponds  to  a  decrease  in  exit  temperature  after 
the  step  input  as  seen  in  Fig.  20.  However,  the  lower  mass  in 
Scenario  3  has  lower  thermal  energy  and  a  greater  drop  in 
temperature  after  the  methanol  step  input.  Furthermore,  the 
effect  of  lag  on  the  decreased  mass  in  Scenario  3  is  twice  that 
of  Scenarios  1  and  2,  indicating  increased  sensitivity  to  a 
delay  in  burner  input  methanol. Fig.  21  is  the  reformer  exit 
hydrogen  response  to  the  methanol  step  input.  There  appears 
to  be  negligible  affect  of  the  lag  condition  on  the  hydrogen 
generated  in  the  reformer,  although  it  appears  that  slightly 
less  hydrogen  is  produced  in  the  lag  condition.  The  drop  in 
hydrogen  corresponds  to  a  drop  in  the  conversion  ratio  since 
the  input  methanol  is  stepped  up  to  the  maximum  level. 
During  this  period  the  thermal  energy  has  yet  to  be  trans¬ 
ferred  fully  from  the  burner  to  the  fuel  in  the  reformer.  For 
the  remainder  of  the  simulation  the  hydrogen  responds  to  the 
maximum  input  methanol  with  first-order  behavior. 

Figs.  22  and  23  are  the  ideal  and  lag  reformer  exit  carbon 
monoxide  responses  to  the  methanol  step  input.  The  lagging 


condition  produces  slightly  less  CO  than  the  ideal  condition 
since  there  are  slightly  less  reactions  to  convert  methanol  to 
CO  and  hydrogen.  This  is  consistent  with  the  plots  in  Figs.  17 
and  21  in  which  there  is  more  unconverted  methanol  and  less 
exit  hydrogen  for  the  lag  condition.  The  CO  drops  with  the 
step  up  in  methanol  flow  and  gradually  reaches  steady  state 
for  the  remainder  of  the  simulation. 

Another  point  with  respect  to  the  CO  response  is  worth 
mentioning.  The  methanol  step  input  is  preceded  by  an  idle 
condition  during  which  the  methanol  input  flow  is  held 
constant  for  around  200  s  and  then  decreases  as  shown  in 
Fig.  24.  This  results  in  an  upward  spike  in  CO,  shown  in 
Fig.  25  for  the  ideal  condition  and  Fig.  26  for  the  lag 
condition,  since  whatever  MeOH  that  exists  in  the  reactor 
converts  to  CO  as  the  reformer  temperature  is  high.  This 
behavior  is  also  seen  in  experiments  conducted  by  [13]. 

The  results  presented  have  illustrated  some  effects  of  lag 
in  burner  input  methanol  on  a  limited  set  of  fuel  processor 
output  parameters,  such  as  unconverted  methanol,  exit 
hydrogen  and  exit  CO.  The  results  also  highlight  conditions 
that  could  exhibit  sensitivity  to  lag  such  as  reduced  thermal 
mass.  The  results  presented  here  have  focused  on  a  meth¬ 
odology  to  illustrate  these  effects  by  using  a  transient  step 
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Fig.  21.  Reformer  exit  H2:  ideal  condition  ( — ),  lagging  condition  (...). 
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Fig.  23.  Reformer  exit  CO — lagging  condition. 
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Fig.  24.  Idle  condition  flowrate  before  step  input:  reformer  ( — ),  burner 
(•••). 
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Fig.  25.  CO  response  to  drop  in  flowrate — ideal  condition. 
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Fig.  26.  CO  response  to  drop  in  flowrate — lagging  condition. 


response.  Preliminary  results  for  a  step  input  show  that  there 
appears  to  be  minimal  effect  of  a  lag  on  the  order  of  2  s. 


3.3.  Dynamic  impact 


Fig.  27.  Transient  step  response  of  the  fuel  processor.  Decreasing  levels  of 
thermal  integration  (from  case  a  to  case  c). 


H2  response  to  a  step  change  in  methanol  flow  from  idle  to 
100%  for  burner/reformer  systems  with  different  thermal 
integration.  Case  (a)  represents  the  case  in  which  Ubmner_ 

wall— side^wall  tffld  ^reformer— wall— side^ reformer  4re  60  ^V/K.  Case 

(b)  represents  a  case  in  which  the  UA  values  are  an  order  of 
magnitude  lower  than  case  (a),  and  case  (c)  UA  values  are  an 
order  of  magnitude  lower  than  case  (b).  It  can  be  seen  that 
the  dynamic  response  time  for  case  (c)  is  around  50  s. 

In  case  (c),  there  is  poor  conversion  of  methanol  to 
hydrogen  during  the  transient  period.  What  actually  happens 
is  that  due  to  poor  thermal  integration  between  the  burner 
and  the  reformer,  the  thermal  energy  from  the  burner  cannot 
be  quickly  transferred  to  the  reformer.  This  results  in  a  drop 
in  the  reformer  catalyst  temperature  during  the  initial  tran¬ 
sient  period  and  this  in  turn  reduces  the  catalyst  activity. 

For  case  (a),  representing  the  highest  relative  level  of 
thermal  integration,  the  result  is  a  quicker  dynamic  response. 
In  contrast  to  cases  (b)  and  (c),  the  response  time  is  1-2  s.  A 
point  to  note  is  that  in  addition  to  thermal  integration, 
thermal  mass  can  play  a  role  as  energy  storage  that  maintains 
a  temperature  at  which  activity  can  occur  and  methanol  can 
be  converted. 

The  nature  of  the  transient  response  also  deserves  some 
explanation.  The  initial  response  from  the  fuel  processor 
seems  instantaneous.  This  is  because  the  reformer  has  a 
finite  volume  and  the  latter  half  of  the  fuel  processor  is  rich 
in  hydrogen.  When  the  methanol/water  mixture  comes  into 
the  reformer,  it  effectively  pushes  the  hydrogen-rich  refor¬ 
mate  out  of  the  fuel  processor  (note  that  the  reformer  is 
operated  in  a  constant  pressure  mode).  However,  this  sce¬ 
nario  cannot  be  sustained  and  the  final  rise  in  the  conversion 
is  due  to  the  fact  that  the  thermal  energy  transfer  approaches 
the  final  steady  state  value.  Note  that  steam  generation  and 
fuel  vaporization  have  an  impact  on  transient  response; 
however,  while  these  processes  are  important  they  have 
not  been  considered  here. 


During  dynamic  operation  of  the  fuel  cell  vehicle,  the 
response  of  the  fuel  processor  becomes  critical  and  is  very 
strongly  influenced  by  the  extent  of  thermal  integration 
between  the  reformer  and  the  burner.  Fig.  27  shows  exit 


3.4.  Burner  emissions 

Two  burner  configurations  have  been  studied:  (a)  adia¬ 
batic  and  (b)  heat  loss  via  thermal  integration  with  the 
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reformer.  For  each  configuration,  temperature  profiles  and 
NOx  formation  are  compared.  The  catalytic  burner  simula¬ 
tion  results  forecast  low  NO*  emissions  from  the  burner  for 
case  (b). 

3.4.1.  Thermal  Integration  effect  on  burner 
temperature 

Fig.  28  shows  the  temperature  profile  of  the  adiabatic  and 
thermally  integrated  cases,  indicating  a  lower  and  more 
uniform  temperature  distribution  in  thermal  integration. 


The  reaction  rate  expression  and  reaction  rate  and  equili¬ 
brium  constants  are  obtained  from  [17].  Deriving  rate 
expressions  for  NO  using  the  reaction  Eqs.  (13)  and  (14) 
and  applying  the  steady  state  assumption  for  N  yields  the 
expression  in  Eq.  (16).  The  subscripts  for  the  reaction  rate 
constants  correspond  to  reaction  Eqs.  (13)  and  (14). 


d[NO]  2fe13[0][N2]{l  -  ( [NO] 2/ Kp^o [N2]  [02] ) } 


d  t 


l  +  (*_13tN0]/*14[02]) 


(16) 


Parameter  values  are  as  follows: 


3.4.2.  NOx  emissions 

It  is  well  accepted  [14-16]  that  there  are  three  types  of 
NOx  formation:  thermal  NOx,  prompt  NOx  and  fuel  NOx. 
Thermal  NOx  is  formed  by  the  oxidation  of  nitrogen  and  is 
nearly  linearly  dependent  on  residence  time  and  rises  expo¬ 
nentially  with  temperature.  Prompt  NOx  occurs  in  reactions 
occurring  under  fuel  rich  conditions.  Fuel  NOx  results  from 
the  oxidation  of  nitrogen-containing  compounds  within  the 
fuel.  The  following  assumptions  were  made  for  the  NOx 
formation  model: 

(1)  Methanol  has  no  nitrogen-containing  compounds  [14] 
and  the  combustion  mixture  is  lean  (stoichiometric 
ratio  =  1.5);  therefore,  only  thermal  NOx,  and  not 
prompt  or  fuel  NOx,  occurs.  The  prediction  of  the  NOx 
emission  characteristics  is  made  using  a  model  in 
which  NOx  formation  occurs  only  in  the  gas. 

(2)  Of  the  dominant  compounds  considered  under  the 
category  “NOx,”  such  as  nitric  oxide  (NO),  nitrogen 
dioxide  (N02),  and  nitrous  oxide  (N20),  NO  comprises 
close  to  95%  of  NOx  emissions  in  combustion 
processes  [17].  Therefore,  only  NO  emissions  are 
modeled.  The  reactions  equations  that  make  up  the 
Zeldovich  mechanism  [18]  for  NO  formation  are  well 
accepted  [14-17]: 

0  +  N2  ^NO  +  N  (13) 

N  +  02  ^  NO  +  O  (14) 

0.5O2  ^  O  (15) 


kl3  =  1.4  x  io14e-315-5/«r  x  10“6, 
&14  =  6.4  x  ioV2615/*r  x  10“6, 
k-u  =  1.6  x  ioV161'5^  x  10"6, 


*13*14  = 


Using  these  equations  and  values  for  Kp  NO  as  provided 
in  [17]  for  N2  +  02  <=$>  2NO,  k_13  can  be  solved.  [O]  in 
Eq.  (15)  is  then  determined  using  the  expression 
tfp,0  =  [0]e(^T)1/2/[02]y2  with  the  assumption  that  [O] 
may  be  set  equal  to  [0]e  (e:  equilibrium)  in  the  hot  products 
for  the  reaction  0.502  <^>  O.  Values  for  Kp  0  are  found  in  the 
literature  [19]. 

Fig.  29  shows  the  equilibrium  NO  concentration  values 
calculated  in  the  model  using  the  above  equations  validated 
against  the  predicted  equilibrium  NO  concentration  values 
as  published  by  [17].  Correlation  between  the  model  values 
and  the  published  predictions  is  greater  than  99%  using  the 
Pearson  correlation  method. 

Fig.  30  shows  the  NO  emissions  for  the  adiabatic  and 
thermally  integrated  cases.  As  Fig.  28  shows,  even  while  gas 
adiabatic  temperatures  are  over  1000  K,  the  emissions  in 
Fig.  30  for  the  adiabatic  case  (a)  with  a  step  input  of  methanol 
are  considerably  lower  than  emissions  formed  at  equilibrium 
at  similar  temperatures.  As  stated  earlier,  thermal  NOx 
is  nearly  linearly  dependent  on  residence  time  and  rises 
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Fig.  28.  Temperature  along  burner:  (a)  adiabatic  case,  (b)  thermally 
integrated  case. 


Temperature,  K 

Fig.  29.  Published  values  and  model  values  for  equilibrium  NO 
concentration:  ( — )  published  values,  ( - )  model  values. 
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Fig.  30.  NO  emissions  for  adiabatic  and  thermally  integrated  cases:  (a) 
adiabatic,  (b)  thermally  integrated. 

exponentially  with  temperature.  A  step  input  flowrate  of 
methanol  implies  a  flowrate  of  nitrogen  and  oxygen  in  the 
combustion  products,  or  a  finite  residence  time;  therefore, 
NO  emissions  are  lower  than  in  equilibrium  in  which  nitro¬ 
gen  and  oxygen  molecules  have  effectively  an  infinite  resi¬ 
dence  time.  This  behavior  is  supported  by  [17],  whose 
studies  indicate  that  it  is  undesirable  to  maintain  the  reacting 
mixtures  at  high  temperatures  for  long  periods  of  time. 

Thermal  integration  case  (b)  appears  to  have  a  significant 
impact  on  the  level  of  NO  emissions  by  lowering  the  burner 
temperature.  Experimental  results  seem  to  support  this 
observation;  for  example,  the  work  by  [5]  using  the  plate 
configuration  for  methane  combustion  resulted  in  NOx 
emissions  below  1  ppm. 

4.  Conclusions 

Section  2  described  the  steps  by  which  a  detailed  model 
was  developed  for  a  thermally  integrated  burner/reformer 
system.  Key  results  for  Section  2  are  listed  below: 

(1)  The  PFR/CSTR  methodology  was  selected  from  the 
literature  to  build  a  numerical  model  that  balanced 
accuracy  with  computation  time. 

(2)  Results  from  simulations  were  analyzed  to  determine 
the  importance  of  added  complexity  in  the  burner 
model  and  revealed  that  the  simpler,  less  computation¬ 
ally  intensive  model,  or  the  pseudo-homogeneous 
model,  could  be  used  in  lieu  of  the  heterogeneous 
model. 

(3)  Key  parameters  were  made  scalable  while  maintaining 
the  validity  of  using  the  original  reaction  kinetics. 

Section  3  described  the  simulation  results  that  show  the 
effects  of  thermal  integration  on:  (a)  reformer  temperature 
and  emissions  especially  when  varying  spatial  parameters 
such  as  catalyst  loading,  (b)  sensitivity  to  changes  in  input 
flow  or  thermal  mass,  (c)  dynamic  response,  and  (d)  burner 
emissions.  Key  results  for  Section  3  are  listed  below: 

(1)  Varied  catalyst  loading  can  play  a  role  in  decreasing 
reformer  temperature  and  CO  emissions. 

(2)  For  the  parameters  used  in  the  sensitivity  analysis,  a 
two  second  lag  in  the  step  input  of  burner  input 


methanol  had  a  more  pronounced  (while  still  minimal 
overall)  effect  on  plate  configurations  that  had  a  higher 
heat  transfer  coefficient  and  lower  thermal  mass  than 
the  base  case. 

(3)  Dynamic  response  can  be  greatly  affected  by  thermal 
integration:  a  thermally  well-integrated  system  can 
have  an  order  of  magnitude  decrease  in  step  response 
times. 

(4)  Reducing  burner  temperature  through  thermal  integra¬ 
tion  can  minimize  burner  emissions  such  as  nitric 
oxide. 

A  final  conclusion  reached  from  this  research  is  that  the 
full  impacts  of  thermal  integration  have  yet  to  be  explored 
and  studied  in  the  context  of  the  fuel  processor  and  the 
overall  indirect  methanol  fuel  cell  system. 

5.  Future  research 

This  work  has  focused  on  results  from  a  fully  warmed  up, 
thermally  integrated  burner/reformer  system  using  standard 
industrial  catalysts.  Further  development  on  the  fuel  pro¬ 
cessor  could  focus  on  the  following  issues: 

•  Cold  start:  warm  up  time,  emissions  and  fuel  consump¬ 
tion. 

•  Different  burner  and  reformer  catalysts. 

•  Other  variations  in  burner  catalyst  loading. 

•  Variations  in  reformer  catalyst  loading. 

•  Control  strategies  that  affect  the  reformer’s  ability  to 
provide  the  necessary  hydrogen  to  the  fuel  cell  stack. 
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Appendix  A.  Validation  of  heterogeneous 
burner  model 

The  heterogeneous  burner  model  was  based  on  the  study 
by  [8],  and  this  study  was  selected  for  the  following 
reasons.  [8]  conducted  experimental  and  numerical  ana¬ 
lysis  to  evaluate  the  performance  of  a  catalytic  combustor 
with  a  methanol  mixture.  Other  studies  discuss  transient 
analyses  with  methane,  propane  and  other  hydrocarbon- 
air  mixtures.  [8]  was  consulted  for  the  present  work 
because  of  the  provision  of  surface  and  gas  reaction  rate 
constants  for  methanol  combustion.  Furthermore,  [8]  pro¬ 
vides  a  detailed  analysis  of  start-up  emissions,  which  is  not 
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considered  in  the  present  work  but  will  be  a  reference  for 
future  work. 

The  prediction  of  the  temperature  profile  is  made  using  a 
heterogeneous  model  in  which  reactions  occur  both  on  the 
surface  and  in  the  gas.  The  following  assumptions  were 
made  in  the  derivation  of  the  energy  and  specie  balance 
equations  for  the  catalytic  combustion  of  methanol: 

(1)  Radiation  within  the  combustor  and  to  the  surroundings 
is  neglected. 

(2)  Pressure  drop  is  neglected. 

(3)  The  heat  loss  from  the  combustor  is  assumed  to  be  only 
to  the  reformer  and  not  to  the  ambient  environment. 

(4)  The  gas  mixture  is  an  ideal  gas. 

(5)  Intermediate  reaction  products,  such  as  formaldehyde 
and  carbon  monoxide,  are  not  considered.  Complete 
conversion  of  methanol  to  carbon  dioxide  and  water  is 
assumed. 


The  modeling  methodology  uses  the  CSTR  approach  and 
energy  and  species  balance  calculations  are  performed  in 
each  CSTR.  The  surface  energy  balance  is  shown  in 
Eq.  (A.l). 


'stored 


=  E 


in 


Peat  ^Pcat  ^catAc  Av 


2  dTcat 
d  t 


ScatAcA(rcat(x+1)  2Tcat(v)  -f-  Tcat(x_1))  T  ScatAcAx 
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(7g  -  Tc at)  +  ]T 


d  rij 


*cat 


d  t 
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+E 


_  prod 
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Finite  difference  formulation  is  used  for  the  heat  conduction 
term,  and  enthalpy  flow  to  and  from  the  surface  and  gas  is 
taken  into  account.  With  the  exception  of  parameters 
described  earlier,  listed  below  are  details  of  the  parameters 
in  Eq.  (A.l). 


•  Scat  is  the  closed  frontal  area  fraction  (i.e.  1  -S,  where  S  is 
the  open  frontal  area  fraction  of  0.64  [8]. 

•  Ac  is  the  reactor  cross-sectional  area,  and  Av  is  the  length 
of  one  CSTR. 

•  In  the  conduction  term,  2  is  the  gamma  alumina  thermal 
conductivity  of  0.94  W/mK  [8],  rcat(x)  is  the  catalyst 
temperature  of  one  CSTR,  Tcat(x+1)  is  the  catalyst  tem¬ 
perature  for  the  CSTR  downstream  of  Tcat(x)  (recall  that 
the  plug  flow  reactor  is  represented  by  a  series  of  10 
CSTRs),  and  Tcat(x_1)  is  the  catalyst  temperature  for  the 
CSTR  upstream  of  rcat(x). 

•  In  the  convection  term,  the  specific  surface  area,  £  is 
2700  m2/m3  [8],  the  heat  transfer  coefficient,  ht  is  94  W/ 
nr  K  [9]  and  Tg  is  the  gas  temperature. 

•  The  third  term  on  the  right  hand  side  of  the  equation 
represents  heat  flow  of  surface  species  arriving  as  reac¬ 
tants  and  leaving  as  products. 


The  gas  energy  balance  is  shown  in  Eq.  (A. 2). 


d7g 
d  t 


-AcAx£ht(Tg  -  :rcat)  - 


prod 


(A. 2) 


The  mass  transfer  coefficient,  PD,  rate  constants  k\  and  k2 
and  reaction  and  rate  equations  are  the  same  as  shown  earlier 
in  Eqs.  (l)-(3).  The  species  balance  equations  are  the  same 
as  shown  earlier  in  Eqs.  (4)-(7).  Only  Eq.  (8)  differs  slightly 
as  shown  in  Eq.  (A. 3):  instead  of  a  single  temperature  Tburner, 
Tg  is  used,  as  the  gas  and  surface  temperatures  are  deter¬ 
mined  separately  in  the  heterogeneous  model. 

=  Flim)  -  y>g  +  r„J  +  (A3) 

Note  that  calculations  are  done  with  moles  and  therefore 
use  the  universal  gas  constant/?;  however,  while  moles  are  not 
conserved,  mass  is  conserved.  This  condition  was  verified 
within  the  model.  The  initial  condition  for  the  surface  and  gas 
temperatures  are  Tcat  =  Tg  =  550  K,  where  the  surface  and 
gas  are  assumed  to  be  fully  warmed  up  at  time  t  =  0.  The 
initial  conditions  for  moles  are  no2  =  0.2ht,  and  n^2  =  0.8^t 
assuming  only  air  exists  within  the  burner  at  t  =  0. 

Fig.  31  shows  a  comparison  of  the  [8]  and  this  paper’s 
model  temperature  profiles  using  parameters  used  in  [8] — 
input  temperature  of  393  K,  relative  air/fuel  ratio  of  3.4,  and 
space  velocities  of  3700,  5700,  8500  and  13000  h_1.  Space 
velocity  (SV)  is  defined  as  the  ratio  of  the  feed  volumetric 
flowrate  at  STP  to  the  reactor  volume.  [8]  provided  no  gas 
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Fig.  31.  Comparison  of  published  and  model  temperature  profiles. 
Relative  air/fuel  ratio  =  3.4,  inlet  T  =  393  K  and  space  velocities  (SV)  are 
as  indicated  [8].  ( — )  [8],  ( - )  this  paper’s  model. 
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temperature  profiles  in  this  particular  study  so  a  comparison 
of  the  catalyst  temperature  is  shown  here. 

Simulation  results  of  the  adiabatic  case  were  compared 
with  the  experimental  temperature  profiles  provided  by  [8] 
with  a  correlation  of  greater  than  98%  using  the  Pearson 
correlation  method. 


Appendix  B.  Discussion  of  UA  parameter 


t 


burner-catalyst 


kef-catalyst 


1/U 


burner-wall -side 


1/U 


re  former-wall -side 


Fig.  33.  Resistances  across  plate. 


The  main  parameter  used  in  the  models  to  thermally 
couple  the  burner  and  reformer  is  UA,  where  U  is  the  heat 
transfer  coefficient  and  A  is  the  area.  Specifically,  UA  is  the 
product  of  the  heat  transfer  coefficient  between  the  burner 
catalyst  and  wall  and  wall  area  (t/bumer-waii-side^waii)  and  the 
product  of  the  heat  transfer  coefficient  between  the  wall  and 
reformer  catalyst  and  wall  area  (t/reformer-waii-side^waii)- 
The  best  way  to  determine  these  UA  values  would  be  from 
microscale  experiments.  It  might  also  be  possible  to  estimate 
UA  from  geometrical  and  thermal  properties  provided  the 
appropriate  physical  properties  can  be  accurately  deter¬ 
mined.  Alternatively,  given  a  value  for  UA,  the  correct  set 
of  thermal  properties  can  lead  to  initial  geometry  estimates. 

Both  methods  require  a  set  of  heat  transfer  equations.  To 
illustrate  the  parameters  used  in  these  equations,  Fig.  32 
shows  one  plate  with  burner  and  reformer  catalyst. 

Note  the  inclusion  of  the  film  resistance  on  the  hot  burner 
side  and  the  film  resistance  on  the  cold  reformer  side.  In  the 
pseudo-homogeneous  reactor,  the  warming  of  gases  and  the 
catalyst  is  taken  into  account,  but  the  gas  is  in  thermal 
equilibrium  with  the  catalyst.  Most  of  the  reactions  take  place 
on  the  surface;  therefore,  most  of  the  heat  is  generated  by  those 
surface  reactions.  However,  if  the  gas  thermal  storage  rate  were 
significant  relative  to  reaction  heat  rate,  the  barrier,  or  film 
resistance  between  the  gas  and  catalyst  would  have  to  be  taken 
into  account,  requiring  a  heterogeneous  reactor.  It  was  found 
that  the  gas  thermal  storage  rate  was  not  significant  relative  to 
the  reaction  heat  rate  in  both  the  reformer  and  burner.  This 
claim  can  be  shown  quantitatively  from  model  results. 

The  thermal  storage  rate  is  calculated  using  (XX^outz- 
hm i)cv.)AT,  where  the  flowrates  are  the  in  and  out  flows  of 


^burner-catalyst 

kvall 

kef-catalyst 


Fig.  32.  Cross-section  of  plate  with  burner  and  reformer  catalyst.  Not  to 
scale;  t  =  thickness. 


Burner  film 


Reformer  film 


Table  1 


Ratios  of  thermal  storage  rate  to  reaction  heat  rate 


Regime 

B  (burner), 

R  (reformer) 

Thermal  storage 
rate  (J/s) 

Reaction  heat 
rate  (J/s) 

Ratio 

Inlet — exit 

B 

1052 

24435 

0.043 

CSTR  1 

R 

529 

14419 

0.037 

each  CSTR  for  the  ith  specie,  ATis  the  temperature  difference 
of  the  input  and  output  flows,  and  cVi  is  calculated  for  each 
specie  at  the  average  temperature.  The  reaction  heat  rate  is 
calculated  by  using  («cH3OHin  -  «ch3ohoui  )  A//s,  where  the  net 
methanol  flowrates  for  each  CSTR  are  multiplied  by  the 
reaction  heat.  The  reaction  heat  for  methanol  oxidation  is 
640,800  J/mol  and  for  methanol  steam  reformation  is 
50,000  J/mol.  Table  1  shows  the  ratios  of  thermal  storage  rate 
to  reaction  heat  rate  for  just  the  first  CSTR,  as  it  is  the  highest 
among  the  series  of  CSTRs  for  both  burner  and  reformer. 

As  the  results  show,  the  ratio  of  thermal  storage  rate  to 
reaction  heat  rate  is  negligible;  therefore,  it  was  decided  to 
ignore  the  film  resistances  and  consider  only  the  burner 
catalyst,  reformer  catalyst  and  wall.  The  resistance  to  heat 
transfer  by  conduction  is  shown  in  Fig.  33.  Note  that  the 
nodes  defining  the  resistance  boundaries  are  placed  at  the 
center  of  each  section.  This  method  was  used  in  order  to  use 
the  average  temperature  in  the  section,  which  was  assumed 
to  be  in  the  center  under  steady  state  conditions. 

B.l.  Equations  used 

Fig.  34  provides  the  context  by  which  the  equations  were 
derived. 

The  parameters  k  and  t  are  the  thermal  conductivities  and 
material  thickness,  respectively.  An  arbitrary,  decreasing 
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Fig.  34.  Diagram  for  derivation  of  heat  transfer  equations. 
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temperature  profile  is  assumed,  represented  by  71,  T2 ,  73, 
T4  and  75.  At  steady  state,  the  heat  flux,  q ",  is  identical  in 
each  section  and  can  be  described  by  the  following  equa¬ 
tions.  For  the  burner-wall-side: 

general  equation  :  q"  =  U AT 


a 


a 


or 


4' 


or 


^burner— catalyst  /  2 

//  ^burner— catalyst/ 2 

^burner— catalyst 


Tl  -  T2 


(B.l) 


 ^wall 


Evall  /  2 

//  Evall/2 

q  ~1 - 

^wall 


(73  -  73) 


T2  -  73 


(B.2) 


Adding  heat  flux  Eqs.  (B.l)  and  (B.2)  together,  U  is  deter¬ 
mined  as  shown  in  Eq.  (B.3). 

// ^burner— catalyst/^  ,  6vall/2\ 

q  (  — - £ - h  — - I  =  7  1  —  7  3 


£/  = 


^burner— catalyst  ^wall  / 

^burner — catalyst  k wall 

(^wall/2) ^burner— catalyst  (^burner— catalyst  /  2) &wa]] 


(B3) 


A  similar  analysis  can  be  made  on  the  reformer- wall  side. 
The  full  UA  equations  are  shown  in  Eqs.  (B.4)  and  (B.5). 


^burner — wall — side  A  wall 

^wall  ^burner— catalyst  (  Eb 


urner- 


-catalyst/  A^bumer— catalyst) 


(^wall/2) ^burner— catalyst  “E  (^burner— catalyst /2)Ewa|] 


(B.4) 


Table  2 

Material  thermal  conductivities 


Plate  section 

Material 

Thermal  conductivity 
(W/m  K) 

Burner 

y  Alumina  [8] 

0.94  [8] 

Wall  [20] 

Stainless  steel  316 

14 

Plain  carbon  steel 

60 

Al  alloy  2024-T6 

177 

Reformer 

Cu/Zn0/A1203  [12] 

15  [21-23] 

If  the  values  of  UA  are  estimated  from  the  bulk  material 
properties  of  the  reformer  and  burner  catalysts  along  with 
three  different  wall  materials  (Table  2),  it  was  found  that 
burner— wall— sideAwaii  is  up  to  an  order  of  magnitude  larger 
than  f/reformer-waii-sideAwalb  a  direct  effect  of  the  much 
smaller  burner  catalyst  mass  (6  kg  burner  versus  30  kg 
reformer).  However,  it  is  unclear  whether  the  properties 
of  the  burner  catalyst  will  be  valid  for  a  coated,  monolith 
geometry  where  the  thermal  path  length  and  cross-sectional 
areas  (when  accounting  for  the  voids)  may  be  different, 
again  indicating  the  need  for  laboratory  studies  for  deter¬ 
mining  thermal  properties  of  these  materials. 

Given  that  published  experimental  data  about  UA  values 
for  integrated  burner/reformer  systems  is  scarce,  as  a  first 
step  Eqs.  (B.4)  and  (B.5)  can  be  used  to  determine  the  initial 
estimates  for  additional  simulations. 
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^reformer — wall — side  A  wan 

_  ^wall  ^reformer— catalyst  (  ^reformer— catalyst  / N ^reformer— catalyst ) 

( Eval  1  / 2) ^reformer —catalyst  ~E  (^reformer— catalyst / 2) kwap 
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The  terms  Vh/Nth  and  VrINtr,  where  V  is  volume  and  N  is 
number  of  plates,  represent  per  plate  areas  of  the  burner  and 
reformer,  respectively  and  are  assumed  to  be  equal  in  the 
plate  configuration.  Another  area  term,  Vw/Ntw  could  be 
used  in  lieu  of  either  of  the  other  two  terms,  as  they  are  all 
equal,  and  given  the  assumption  that  the  heat  transfer  area  is 
the  same  as  the  physical  area. 

In  the  initial  part  of  the  analysis,  the  values  for  UA  were 
selected  using  the  work  done  earlier  by  [7].  The  study  by  [7] 
on  a  shell  and  tube  reformer/burner  configuration  illustrated 
the  need  for  substantially  higher  heat  transfer  coefficients  to 
achieve  exit  H2  response,  from  10  to  90%,  of  the  order  of  2  s, 
which  is  the  response  time  required  by  the  fuel  cell  system 
specifications  used  [11].  Therefore,  the  initial  studies  of 
the  thermally  integrated  reformer/burner  system  were  based 
on  an  assumption  of  setting  both  f/bUmer-waii-sideAWaii  and 
^former— wall— sideAwaii  equal  to  60  W/K  per  plate  (for  100 
plates),  which  is  an  order  of  magnitude  higher  than  the 
values  studied  by  [7]. 
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